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Introduction 


A.  Origin  of  Problem 


Suppose  there  come  into  your  possession  one  or  more 
"black  boxes"  and  the  information  that  these  are  linear.  Why 
would  you  want  to  determine  their  impulse  responses,  and  how 
would  you  go  about  it? 

A fairly  literal  case  would  be  that  you  had  purchased 
the  systems  "impulsively"  at  a surplus  shop,  and  discovered 
that  the  innards  were  potted  so  as  to  make  direct  inspection 
impossible.  It  would  then  be  desirable  to  measure  the  impulse 
response  or  the  transfer  function  (its  Fourier  transform)  of 
each  system  to  determine  how  it  could  be  utilized. 

In  a ballistic  missile  defense  (BMD)  context  the  linear 
systems  in  question  may  model  how  projectiles  of  interest  to  the 
defense,  re-entry  vehicles  (RV's),  decoys,  chaff  particles,  etc. 
respond  in  time  to  the  voltage  waveform  transmitted  by  a radar. 
One  of  the  defense  functions  is  discrimination  - the  task  of 
determining  which  targets  must  be  acted  against  and  which  are 
there  to  confuse  or  saturate  the  defense.  One  of  the  methods 
for  doing  this  is  to  utilize  differences  in  the  size  or  shape 
of  the  targets.  Direct  inspection  is  again  impossible,  so  the 
received  radar  waveform  is  used.  Sizing  a target  along  the 
radar  line  of  sight  amounts  to  explicitly  or  implicitly 
determining  and  then  further  operating  on  the  impulse  response 
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of  a linear  system.  If  the  processing  technique  extracts  an 
estimate  of  the  target  extent  along  the  l^ne  of  sight  it  is 


referred  to  as  "length  measurement."  If  the  processing  makes 
more  implicit  use  of  the  impulse  response  it  is  a particular 
application  of  a technique  called  "pattern  recognition."  If 
sizing  is  also  carried  out  in  the  cross-range  direction  using 
Doppler  processing  the  technique  is  referred  to  as  "imaging." 

We  will  restrict  ourselves  in  what  follows  to  the  one-dimensional 
case . 

B.  Deterministic  Signal  Approach 

Having  decided  that  you  want  to  determine  an  impulse 
response,  how  do  you  proceed?  The  basic  procedure  is  to  drive 
the  system  with  a test  signal,  observe  the  response  and  process 
it,  the  processing  depending  on  the  test  signal.  The  most 
obvious  approach  is  to  "pulse"  the  system.  A deterministic 
pulse  narrow  with  respect  to  the  system  response  time  (wideband 
with  respect  to  the  system  bandwidth)  is  applied,  and  the 
resulting  output  signal  is  taken  as  an  estimate  of  the  impulse 
response,  albeit  somewhat  smeared  and  distorted.  Choosing  an 
appropriate  pulse  bandwidth  requires  prior  knowledge  about  the 
system  or  a trial  and  error  procedure. 

In  principle,  the  effects  of  finite  pulse  bandwidth  and 
detailed  pulse-shape  can  be  removed  either  in  the  time  or 
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frequency  domains,  yielding  the  "true"  impulse  response.  In 
practice  the  degree  to  which  this  can  be  done  is  limited  by 
interfering  noise,  and  while  the  effect  of  interfering  noise  is 
the  main  subject  of  this  report,  it  is  assumed  that  this  pro- 
cessing refinement  is  not  involved. 

In  the  frequency  domain  the  unknown  transfer  function 
could  be  traced  out  by  applying  a sinusoid  to  the  system  and 
varying  its  frequency.  The  use  of  a pulse  waveform  can  be 
thought  of  as  a way  of  simultaneously  applying  frequency 
components  to  the  system.  This  report  will  deal  with  the  time 
domain  version  of  the  problem. 

What  if  we  allow  the  simple  pulse  to  be  less  simple? 

That  is,  what  if  we  employ  a signal  having  the  desired  bandwidth 
but  a time-bandwidth  product  greater  than  unity?  The  basic 
reasons  for  employing  such  a waveform  in  radar  practice  are: 

(i)  the  presence  of  interfering  noise  originating  within  the 
radar  receiver  and  possibly  also  from  electronic  countermeasures 
(ECM) , external  jamming  sources  employed  by  the  offense  to  screen 
targets  from  the  radar  (ii)  that  the  achievable  signal-to-noise 
ratio  or  radar  sensitivity  depends  on  how  much  energy  can  be 
packed  into  the  radar  signal  (iii)  and  that  radar  transmitters 
operate  under  a peak  power  level  constraint. 

Thus  a modulation  such  as  linear  frequency  modulation 
(LFM) , or  binary  phase  shift  coding  is  utilized  to  increase  the 
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signal  duration  and  hence  energy  while  maintaining  the  bandwidth 
needed  for  resolution.  The  receiver  processing  must  now  change 
from  a simple  observation  of  the  returned  pulse  shape  for  two 
reasons . 

Firstly,  this  pulse-shape  no  longer  looks  almost  like 
the  test  impulse  response,  due  to  the  longer  signal  duration. 

To  correct  this  the  signal  bandwidth  must  be  utilized  to  achieve 
its  inherent  time  resolution.  The  conventional  way  of  doing  this 
is  to  convert  the  radar  signal  into  a signal  having  a unity 
time-bandwidth  product  by  matched  filtering,  i.e.,  by  convolving 
the  received  signal  with  a time-reversed  replica  of  the  trans- 
mitted signal.  It  is  well  known  that  matched  filtering  generates 
the  autocorrelation  function  of  the  signal  in  question,  which  in 
fact  has  a central  "spike"  whose  time  duration  is  the  reciprocal 
of  the  signal  bandwidth.  The  overall  effect  is  that  of  having 
transmitted  a "narrow  pulse,"  the  autocorrelation  function  of 
the  signal,  and  taken  the  response  of  the  unknown  system  to  this 
waveform  as  an  estimate  of  its  impulse  response. 

Secondly,  the  receiver  processing  must  take  into  account 
the  interfering  noise  which  dictated  the  use  of  a time  extended 
modulated  signal  from  an  energy  standpoint.  Basically  the 
receiver  will  reject  noise  outside  the  signal  bandwidth,  and 
secondarily  it  will  employ  an  optimized  passband  shape  to 
maximize  the  signal- to-noise  ratio.  In  other  words,  it  will 
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employ  the  same  matched  filter  that  was  required  to  achieve  the 
desired  resolution  from  the  radar  signal. 

The  above  paragraphs  have  described  how  a conventional 
pulse  compression  radar  would  be  used  to  measure  an  impulse 
response.  In  the  analysis  to  be  described  in  Section  II,  a 
scalar  or  baseband  (rather  than  complex  or  bandpass)  model  is 
used  for  the  time  domain  signals  and  systems.  Thus  in  the 
pulsed  case  the  details  of  a phase  modulated  pulse  compression 
radar  are  not  modeled.  Rather  the  transmitted  signal  is 
thought  of  as  a "narrow"  pulse  whose  time-bandwidth  product  is 
approximately  unity,  and  the  receiver  as  employing  a filter 
matched  to  that  signal.  This  simplification  in  no  way  compromises 
the  applicability  of  the  results  to  a "real"  radar  system. 

C.  Random  Signal  Approach 

Let  us  now  shift  to  a consideration  of  random  signals 
and  their  possible  application  to  impulse  response  determina*- . ^n. 
As  we  have  noted,  the  use  of  a pulse  signal  can  be  thought  of 
as  a way  of  simultaneously  applying  many  frequency  components  to 
trace  out  the  unknown  transfer  function.  Can  a random  process 
be  used  to  perform  this  function,  and  if  so  why  would  it  be 
desirable  to  use  one?  The  familiar  stationary  "whitish" 

Gaussian  process  so  often  employed  in  radar  analysis  can  be 
assigned  a power  spectrum  which  extends  considerably  beyond  the 
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passband  of  the  unknown  transfer  function  and  is  relatively  flat 
within  the  passband.  Thus  in  an  ensemble  sense  the  random 
process  can  apply  frequency  components  to  trace  out  the  transfer 
function . 

The  power  spectrum  of  the  output  process  is  a tracing 
of  the  magnitude  squared  of  the  unknown  transfer  function.  This 
would  be  sufficient  to  determine  a real,  even  impulse  response, 
which  has  a real  (and  even)  Fourier  transform.  A more  arbitrary 
real  impulse  response  with  both  odd  and  even  components  will 
correspond  to  a complex  valued  transfer  function  whose  magnitude 
and  phase  must  be  traced  out.  The  appropriate  tracing  is  the 
cross  power  spectral  density  between  input  and  output.  Since 
this  function  is  the  product  of  the  input  power  spectrum  and  the 
unknown  transfer  function,  the  use  of  a random  signal  with  a 
whitish  power  spectrum  will  produce  a cross-spectrum  which 
essentially  replicates  the  transfer  function  in  magnitude  and 
phase . 

The  corresponding  relationship  in  the  time  domain  is 
that  the  cross-correlation  function  between  input  and  output  of 
the  impulse  response  is  equal  to  the  input  autocorrelation 
function  convolved  with  the  impulse  response.  This  provides 
the  basis  ^ ^ ensemble  sense  for  determining  the  impulse 
response  using  a random  signal.  If  the  autocorrelation  function 
of  the  random  signal  is  narrow  relative  to  the  impulse  response 
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"duration"  the  cross-correlation  function  is  a slightly  smeared 
tracing  of  the  impulse  response. 

In  practice,  what  would  be  available  for  impulse 
response  determination  is  a finite  time  segment  of  a particular 
sample  function  of  the  random  process.  Operationally,  an 
"ergodic"  assumption  would  be  made,  and  the  time  cross-correlation 
function  between  the  available  segments  of  the  input  and  output 
sample  functions  would  be  treated  as  the  impulse  response 
estimate . 

D.  Relationship  of  Approaches 

We  can  now  note  a striking  similarity  with  the  use  of 
the  deterministic  signal.  In  either  case  the  received  signal 
is  cross-correlated  with  a replica  of  the  transmitted  signal. 

The  resulting  waveform  is  then  the  desired  impulse  response 
smeared  by  the  auto-correlation  function  of  the  transmitted 
signal.  So  what  then  is  the  difference  between  using  determin- 
istic and  random  signals?  In  our  view,  there  really  isn't  any, 
and  there  are  several  ways  of  looking  at  it. 

Suppose  first  that  a given  segment  of  random  process 
sample  function,  once  generated,  is  utilized  over  and  over 
again  every  time  the  radar  attempts  an  impulse  response 
determination.  In  this  case  the  fact  that  the  radar  signal  was 
originally  generated  by  a random  process  becomes  irrelevant. 
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It  is  now  a deterministic  signal.  In  practice  its  .;ime  duration 
would  be  long  relative  to  its  correlation  time,  and  the  operation 


of  correlating  it  against  a replica  of  itself  would  "pulse 
compress"  it  to  a narrow  spike  whose  width  was  the  correlation 
time,  riding  above  a sidelobe  structure.  The  whole  operation 
would  closely  parallel  how  things  would  go  with  a signal  that 
had  been  obtained  by  deterministic  means.  A binary  phase-coded 
pulse  compression  signal  for  example  is  really  just  a stylized 
version  of  a sample  function  of  a continuous  random  process  in 
which  the  process  can  change  states  abruptly  after  every 
"correlation  interval"  rather  than  continuously.  The  fact  that 
such  signals  are  often  referred  to  as  "pseudo-random"  or  "noise- 
like" reinforces  the  comparison. 

The  other  possibility  is  that  a segment  of  a new  sample 
function  is  utilized  each  time  the  radar  makes  an  impulse 
response  determination.  There  would  now  be  some  variability 
due  to  the  signal  itself,  but  not  a variability  that  seems 
significant.  Stated  differently,  "most"  sample  function 
segments  are  equivalent  if  their  time  durations  are  long 
relative  to  the  correlation  time.  That  is,  most  sample  functions, 
when  correlated  against  themselves,  will  produce  an  essentially 
identical  correlation  function.  There  will  be  minute  differences 
in  the  exact  shape  of  the  central  spike  and  the  details  of  the 
sidelobes,  but  in  most  cases  the  resulting  impulse  response 

8 


w 


r ’ 

> 

estimate  will  be  essentially  unaffected. 

If  this  line  of  heuristic  reasoning  seems  at  all 

convincing,  one  must  ask  why  the  notion  has  persisted  that 

the  use  of  a random  signal  has  an  advantage  in  impulse  response 

determination.  The  answer  is  basically  as  follows.  In  the 

random  signal  case  the  operation  was  explicitly  recognized  as 

a cross-correlation  operation.  This  led  to  the  idea  that  any 

interference  which  was  "uncorrelated"  with  the  random  signal 

would  be  suppressed  when  the  signal  and  interference  were 

( *) 

I cross-correlated.  The  textbook  by  Y.  W.  Lee  popularized 

j.  this  notion  as  far  back  as  1960.  There  are  several  problems 

[ with  his  discussion  of  this  topic.  Foremost  is  the  fact  that 

i he  apparently  did  not  recognize  that  the  "conventional"  or 

deterministic  signal  method  was  really  just  as  much  of  a cross- 
correlation method  as  the  "novel"  or  random  signal  method, 
s 

I Secondly,  he  did  not  really  deal  with  the  fact  that  an  uncorre- 

I 

> lated  disturbance  will  not  completely  vanish  when  subjected  to 

a finite-time  cross-correlation.  The  degree  to  which  the 
interference  would  be  suppressed  was  not  examined  and  a 
comparison  with  the  "old"  method  was  not  made. 

* 

Y.  W.  Lee,  Statistical  Theory  of  Communication 
j (John  Wiley  and  Sons,  New  York,  1960 ) . 


The  purpose  of  this  introduction  has  been  to  engender 
an  expectation  that  if  the  two  methods  were  compared  on  the 
basis  of  equal  utilization  of  resources  such  as  energy 
and  bandwidth,  they  would  have  essentially  the  same  performance. 
A comparison  has  been  carried  out  analytically,  described  in  the 
following  sections,  which  bears  out  this  expectation. 


Before  proceeding  with  an  analysis  of  each  method  and 
a comparison  between  them,  we  must  establish  a performance 
measure  which  will  be  the  basis  for  comparison.  Using  either 
method,  the  impulse  response  estimate  at  any  point  on  the  true 
impulse  response  will  be  a random  variable.  In  the  deterministic 
signal  method  this  is  due  to  the  interfering  noise  process,  and 
in  the  random  signal  method  it  is  due  both  to  the  interference 
and  to  the  randomness  of  the  signal  itself.  Thus  the  most 
general  description  of  the  estimators  would  be  their  multi- 
dimensional probability  density  functions  at  n points  along  the 
impulse  response.  If  the  ultimate  goal  is  to  describe  how  some 
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discriminant  which  operates  on  the  estimated  impulse  response 
performs,  the  number  and  location  of  the  points  at  which  the  nth 
order  density  functions  are  evaluated  would  depend  on  the 
particular  discriminant.  For  a length  measurement,  a few  points 
near  the  front  and  back  of  each  vehicle  class  might  suffice  if 
the  signal/interference  ratio  were  above  a certain  value.  For 
a pattern  recognition  approach,  points  spread  across  the  full 
extent  of  each  impulse  response  might  be  necessary. 

We  do  not  want  to  consider  specific  discriminants  in 
this  study.  To  avoid  doing  so,  we  choose  as  a performance 
measure  the  normalized  variance  of  the  impulse  response 
estimates  at  any  point  along  the  response.  The  notion  is  that 
any  discriminant  will  benefit  if  this  measure  is  small.  We 
are  saying  that  we  want  "most"  of  the  random  estimates  produced 
by  a given  estimator  to  lie  in  a "narrow"  (in  a percentage 
error  sense)  band  about  the  true  impulse  response  pattern.  We 
don't  particularly  care  whether  an  estimated  pattern  looks  as 
shown  in  Figure  la,  wherein  the  estimation  errors  quickly 
decorrelate  across  the  pattern,  or  as  shown  in  Figure  lb, 
wherein  the  errors  are  highly  correlated  across  the  pattern. 

Thus  our  performance  measure  is  a first  order  statistic. 
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Fig.  1.  Extremes  of  estimation  error  time  behavior 


2 . Model  Considerations 

In  what  follows , expressions  are  simplified  if  we 
calculate  the  normalized  variance  at  certain  "typical"  points 
on  the  pattern,  although  the  calculation  can  be  made  at  any 
desired  point.  Simplification  is  also  gained  if  the  impulse 
response  is  assumed  to  have  even  symmetry.  Again  the  restriction 
is  made  just  for  convenience.  The  analysis  has  been  carried 
out  for  continuous  impulse  responses  described  by  any  one  of  a 
number  of  unimodal  shapes.  These  include  Gaussian,  sin  x/x, 
cusped  exponential,  and  rectangular.  The  same  shapes  were 
used  to  describe  the  deterministic  measurement  signal  and  the 
auto-correlation  function  of  the  random  measurement  signal. 

The  functional  form  of  the  results  does  not  depend  on  the 


shape  chosen  in  the  operating  regime  of  interest  (see  below) . 

Only  form  factors  such  as  tt,/T,  etc.  change.  The  specific  j 

results  in  this  report  will  correspond  to  the  Gaussian  shape  i 


assumption,  and  the  typical  point  at  which  they  are  evaluated 


is  the  peak  of  the  impulse  response  so  shaped. 

In  addition  to  the  basic  results  for  an  impulse 
response  described  by  a smooth  curve,  some  results  are 
presented  for  an  impulse  response  described  by  either  one  or  a 
pair  of  impulses.  The  single  impulse  case  corresponds  to  a 
point  target  model,  and  the  results  provide  some  insight  in  the 
continuous  case.  The  impulse  pair  model,  for  which  results 
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are  similar,  corresponds  to  a situation  in  which  resolution  can 
be  achieved  between  target  features  in  an  operational  sense,  but 
the  bandwidth  of  the  received  signal  is  still  determined  by  the 
transmitted  signal  bandwidth,  not  the  target  bandwidth.  Such  a 
model  applies  to  a vehicle  which  generates  radar  returns 
primarily  from  a tip  and  a base  such  that  the  "rise  time"  of 
the  signal  is  unaffected  by  the  target  at  any  signal  bandwidth 
being  considered.  The  typical  point  selected  in  the  point  target 
case  is  the  location  of  the  single  impulse.  In  the  impulse  pair 
model  two  typical  points  are  selected,  one  at  an  impulse 
location  and  one  midway  between  impulses.  Figure  2 illustrates 
the  three  cases. 

3 . Operating  Regime 

The  results  which  follow  below  correspond  to  a certain 
operating  regime  for  some  of  the  basic  time  and  bandwidth 
parameters  involved.  This  regime  is  defined  by  the  requirements 
that: 

(i)  the  resolution  of  the  estimate  be  somewhat  finer 
than  the  time  scale  of  the  impulse  response 

(ii)  the  estimate  be  generated  over  a significant 
portion  of  the  impulse  response  extent. 

In  our  basic  case  of  a continuous  impulse  response  the 
first  requirement  amounts  to  saying  that  the  signal  bandwidth 
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TYPICAL  POINTS 


Fig.  2.  Basic  impulse  response  forms. 
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must  be  somewhat  larger  than  the  bandwidth  of  the  system  being 
measured.  If  we  denote  suitably  defined  signal  and  system 
bandwidths  by  W and  Vl^  respectively  (subscript  ”h"  for  impulse 
response  h(t)),  then  we  require  that  W>>W^,  where  a reasonable 
operational  definition  of  >>  is  >10  W^.  Rather  than  describing 
the  system  by  a bandwidth  Wj^,  we  can  describe  it  by  a suitably 
defined  impulse  response  extent  Tj^,  where  1.  The 

requirement  then  becomes  WT|^>>1.  in  our  impulse  pair  model, 
the  system  bandwidth  Wj^  is  infinite  but  is  still  meaningful 
and  the  requirement  to  cleanly  resolve  the  pair  of  impulses  is 
again  expressible  as  WT^  >>  1. 

The  second  requirement  concerns  the  amount  of  data  that 
must  be  processed  and  is  most  easily  described  from  the  viewpoint 
of  the  random  signal  method.  Imagine  the  time  origin  to  be 
centered  under  the  impulse  response-  The  estimate  of  h(o)  is 
then  the  cross-correlation  between  the  transmitted  and  received 
waveforms  at  zero  shift  (propagation  delays  having  been  removed) . 
To  generate  the  estimate  a distance  T^  from  the  origin  requires 
shifting  the  data  by  Th-  If  the  signal  is  transmitted  and  the 
data  sample  collected  over  +T,  estimates  can  be  generated  only 
out  to  +T  along  the  impulse  response.  Thus  at  the  very  least, 
T=Th.  Further,  unless  the  time  window  is  somewhat  longer  than 
this  the  data  falloff  will  be  excessive  and  the  window  function 
will  unduly  affect  the  estimate  at  shifts  approaching  Tj^.  It 
is  clearly  desirable  to  require  that  T >>  or  Wj^T  >>  1. 
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It  then  follows  from  the  first  requirement  that  WT>>1  will  also 
characterize  the  operating  regime.  The  estimation  accuracy 
expressions  that  result  from  an  analysis  taken  to  these  limits 
on  the  parameters  will  show  that  desirable  parameter  settings 
from  an  accuracy  viewpoint  are  consistent  with  the  regime 
described  here. 

B.  Random  Signal  Method 

The  random  signal  method  can  be  described  up  to  a point 
without  specific  reference  to  the  form  of  the  impulse  response. 
The  transmitted  signal  is  a finite  time  segment  of  a sample 
function  x(t)  of  a real  zero  mean  Gaussian  random  process.  The 
signal  exists  for  |t|  £ T,  and  a different  sample  function  is 
used  each  time  the  radar  makes  a transmission.  The  power 
spectrum  of  x(t)  is  assumed  to  be 

S^(f)  = exp[-(f^/2w2)  ] (1) 

with  peak  spectral  height  The  bandwidth  parameter  W has 

been  chosen  as  the  "standard  deviation"  of  the  Gaussian  shape 
of  the  spectrum.  The  corresponding  autocorrelation  function  of 
x(t)  (Fourier  transform  of  Sjj(f))  is 

E [x  (t)  X (t+x)  ] = ^xx^^^  ~ exp[-(2iT^W^T^)  ] (2) 
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where  E[*l  denotes  expectation.  For  large  W,  x(t)  approaches 
white  noise  of  spectral  height  and 

Rxx^"^^  ^ NxUo(t)  (3) 

where  u^(*)  is  the  unit  impulse.  The  average  power  in  the 
process  x(t)  is 


R^^(o)  = /m  WNx  (4) 

so  that  the  average  energy  in  a segment  of  x(t)  2T  seconds 
long  is 

Ex  = 2 /Hr  WTNx  (5) 

a relation  that  will  be  used  later 

The  reflected  signal  y(t)  is  the  segment  of  x(t) 
convolved  with  the  unknown  impulse  response  h(t) : 

T 

y(t)  = y*  x(C)  h(t-^)  dC  (6) 

-T 

The  received  signal  is 

z(t)  = y(t)  + n(t)  (7) 
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where  n(t)  is  the  interfering  noise  process,  assumed  to  have 
the  same  statistics  as  x(t),  to  have  spectral  height  Nj^  rather 
than  Nxf  and  to  be  statistically  independent  of  x(t)  so  that 


Efx'^Ct)  n^(t+T)]  = E[x^(t)]  E[n^(t+T)],  any  k,j,T  (8) 

The  estimator  of  h(t),  found  by  cross-correlating 
z(t)  with  the  transmitted  segment  of  x(t)  is 

T 

h(T)  = — f x(t)  z(t+T)  dt  (9) 

-T 

where  propagation  delay  to  and  from  the  target  has  been  taken  as 
zero.  To  calculate  the  mean  and  standard  deviation  of  h(T)  we 
must  make  some  assumptions  about  h(T). 

1 . Continuous  Impulse  Response 

Assume  that  h(T)  is  a smooth  curve  qualitatively 
described  by  an  effective  time  duration  T^^  and  corresponding 
bandwidth  where 

W^Tj^  ~ 1 ( 10 ) 
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Then  using  Equations  (6)  through  (9)  one  can  show  th^t 


E[h(T)  ] = 


2TN 


X 

I 


T-5 

/ 


R (v)h(v+T)  dv  (11) 


C=-T  V=(-T-5) 


If  we  then  assume  that  is  narrow  with  respect  to  h,  i.e., 


that 


= W/Wj^  >>  1 


then  (3)  can  be  used  to  obtain 

E [h(T)  ] h(T) 

so  that  h is  an  unbiased  estimator  of  h. 

The  variance  of  h is 

Var[h(T)]  = E[(h(T)  - hCt))^] 


(12) 


(13) 


(14) 


One  proceeds  by  plugging  (9)  and  (11)  into  (14).  Appendix  A 
describes  what  comes  next.  It  involves  such  things  as  recognizing 
that  a product  of  two  multiple  integrals  can  be  written  as  a 
higher  order  integral  by  introducing  additional  dummy  variables, 
the  use  of  (8) , and  the  use  of  the  "moment  theorem  for  Gaussian 
variables,"  which  allows  one  to  express  the  expected  value  of 
a product  of  three  or  more  variables  as  a sum  of  expectations 
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of  products  of  two  variables  at  a time.  The  result  is  that  the 
desired  variance  can  be  written  as: 


Var[h(T)]  = Var  [h(T)]  + Var  [h(T)] 
A n 


(15) 


The  Varj^  term  depends  only  on  hCx)  . It 

represents  the  variability  in  the  estimate  due  to  the  randomness 
of  the  transmitted  signal,  and  persists  even  if  there  is  no 
interference  whatsoever.  The  Var^  term  depends  only  on 
and  the  autocorrelations  of  the  signal  and  the  interfer- 

ence. It  corresponds  to  the  fact  that  when  the  interference 
is  cross-correlated  against  the  transmitted  process  for  a finite 
time  the  result  has  an  expected  value  of  zero,  but  a non-zero 
variance.  The  two  terms  are  described  in  more  detail  below. 

a . Vary  Term 

If  we  again  set  WTj^>>l,  assume  for  convenience  that  h(T) 
is  even,  and  evaluate  the  variance  at  the  typical  point  t=0,  we 
obtain 


T-f, 


Var  [h(o)  ] = — j- 
2T 


f f h‘(^))  dv  d^ 

4=-T  V=-T-5 


(16) 


Next  assume  that  h(T)  is  also  Gaussian  shaped  and  write  it  as 


h(T)  = h^  exp[- (27rWj^T)  ] 


(17) 
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With  this  convention,  the  power  transfer  function  of  the  target 
(magnitude  squared  of  Fourier  transform  of  h(T))  is 


lH(f)l^  = exp[-(f^/2W^^) ] , 


A comparison  with  (1)  shows  that  the  bandwidth  parameters  W and 
Wj^  have  been  defined  consistently.  That  is,  the  random  pro- 
cesses x(t)  and  n(t)  can  be  thought  of  as  white  noise  passed 

2 

through  a filter  with  a Gaussian  shaped  passband  I 

whose  "variance  parameter"  is  W^,  while  the  target  is  a filter 

2 

with  Gaussian  shaped  passband  |H(f)  | whose  "variance  parameter" 
2 

IS  . The  generic  parameter  Th  used  above  can  be  thought  of 
as  2-  Using  (17)  in  (16)  and  setting  Wj^T  >>  1 results  in 

a normalized  variance  given  by 


Var^ [h (o) ] 

1 

1 

h^(o) 

2/77? 

W.  T 
h 

where  the  first  factor  would  vary  if  the  definition  of  bandwidth 
or  the  shape  assumed  for  h{T.)  changed,  and  the  second  factor 
expresses  the  fundamental  dependence  involved.  Equation  (19) 
has  an  interesting  interpretation  that  will  be  deferred  to 
Section  B-2  below. 


b.  Var  Term 
n 


This  term  results  when  the  procedure  described  just 
before  Equation  (15)  is  carried  out.  It  corresponds  to  the 
expectation 


If  (2)  is  then  applied  and  the  result  evaluated  for  WT  >>  1 we 
arrive  at  the  normalized  variance 


Var^ [h{o) ] 
h^  (o) 


/tT  w 
2 T 


N 


(22) 


where  we  have  let  hQ=l*  This  entails  no  loss  of  generality  if 
we  think  of  as  representing  the  transmitted  power  level 
modified  by  radar  range  equation  factors,  one  of  which  is 
target  cross-section  level,  i.e.,  h^.  Finally,  if  we  apply 
(5) , the  result  can  be  written 


/\ 

Var^ [h (o) ] 

/2  TT 

h^  (o) 

Ev/N 
^ n 

(231 


These  results  can  be  interpreted  readily  for  either  form  (22)  or 
(23).  To  understand  (22)  one  can  go  back  to  (20)  and  think  of 
the  integral  as  a sum  of  independent  random  variables,  4TW  in 
number,  spaced  by  in  time,  each  being  the  product  of  a 

sample  of  x,  a sample  of  n,  and  the  spacing  factor  1/2W. 

Squaring  the  sum,  taking  the  expectation,  and  using  (8)  results 
in  a new  summation  in  which  each  term  contains  the  expected  value 
of  the  square  of  a sample  of  x times  the  square  of  a sample  of  n. 
Applying  (4) then  results  in 
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which  agrees  with  (22)  to  within  a factor  of  /ii,  the  difference 

corresponding  to  the  fact  that  the  integral  form  is  evaluated 

for  the  actual  Gaussian  shape  of  the  correlation  functions  while 

the  summation  form  uses  a rectangle  approximation.  If  we  then 

2 

write  the  factor  W/T  as  W^/TW  we  see  that  Var^  is  directly 

2 

proportional  to  W , inversely  proportional  to  TW,  and  inversely 

proportional  to  These  factors  are  respectively  the 

bandwidth  dependence  of  the  variance  of  each  term  in  the 

original  sum,  the  dependence  of  the  number  of  independent 

samples  on  time  and  bandwidth,  and  the  power  ratio  of  the 

"signal"  and  interference  samples.  These  identifications 

constitute  our  interpretation  of  (22). 

Regarding  form  (23),  the  inverse  dependence  on  a "signal 

energy  to  noise  power  density"  ratio  is  not  surprising  since  we 

are  dealing  with  the  normalized  variance  of  a signature  amplitude. 

2 

The  increase  of  variance  as  W results  from  the  fact  that  we  are 
operating  with  W>>W^  to  achieve  resolution.  One  factor  of  W 
corresponds  to  a loss  of  energy  outside  the  nominal  target 
bandwidth  W^.  The  other  factor  of  W corresponds  to  a mismatch 
loss  at  the  receiver  because  the  receiver  must  operate  with 
bandwidth  W to  achieve  the  resolution  while  the  power  bandwidth 
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of  the  received  signal  is  essentially  W^. 


This  explanation  is 


repeated  more  thoroughly  in  Section  II  C-1  below. 


2 . Impulsive  Impulse  Response 

In  this  section  we  no  longer  think  of  h{t)  as  being  a 
smooth  curve  but  rather  assume  that  it  is  composed  of  either  one 
or  two  impulses.  As  discussed  in  Section  II  A-2,  the  single 
impulse  represents  an  idealized  point  target,  and  provides  a 
comparison  with  the  continuous  case.  The  two  impulse  model  is 
quite  similar,  and  represents  a target  that  returns  from  a "tip" 
and  a "base"  with  essentially  no  broadening  of  either  return  at 
any  operating  bandwidth  being  considered.  The  two  models  are 
shown  in  Figure  3a.  For  simplicity  the  two  impulses  are 
assumed  to  have  equal  areas  and  to  sit  symmetrically  about  the 
time  origin. 

Equations  (1)  through  (9)  and  (11)  still  apply. 

Studying  (11)  indicates  that  in  the  point  target  case  the 
expected  value  of  the  estimate  follows  the  shape  Rxx^^^'  '^hich 
is  a "narrow"  continuous  pulse  with  width  inversely  proportional 
to  W rather  than  a spike.  In  the  two  impulse  case,  E[h(T)] 
looks  essentially  like  a pair  of  offset  replicas  of  ^xx ^ ^ ^ • 

If  1/W  is  small  compared  to  the  impulse  spacing,  i.e.,  if  WTj^>>l, 
the  two  peaks  are  well  resolved  and  E[h(T)]  looks  like  a pair 
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of  "narrow"  continuous  peaks.  Figure  3b  shows  these  curves. 

If  we  proceed  from  (11)  using  (2)  and  assume  that  WTj^  >>  1 
and  WT  >>1,  we  obtain 

E[h(T)]  = /Tn  WhQ,  T at  an  impulse  location 

(25) 

/s 

E[h(T)]  =0,  T "far"  from  an  impulse  location 

where  hQ  is  the  area  of  each  impulse.  The  first  of  these 
values  will  be  used  to  normalize  the  variance  terms  to  follow. 

a.  Varjj  Term 

i 

Equation  (15)  still  applies,  stating  that  the  variance 

of  the  estimate  is  the  sum  of  a Varv  and  a Var  term.  However, 

A n 

the  Varj^  term  is  no  longer  given  by  (16).  The  problem  must  be 
reworked  starting  with  (11)  and  (14)  as  in  Appendix  A.  If  this 
is  done  utilizing  (2)  and  the  model  of  Figure  3a,  assuming  that 
WTh  and  WT  are  >>  1,  and  letting  h^  be  unity,  we  obtain  in  the 
point  target  case 


Varj^[h(T)  ] 

1 

at  impulse  location 

(E[h(o)  ])^ 

2 /ttWT 

1 

4/TrWT 

"far"  from  impulse 
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In  the  two  impulse  case  we  obtain 


i 

t 


f 


I 

I 


I 


i 


Var^[h(T) ] 3 

„ = ;;;; — at  impulse  locations 

(E  [h(Ti)  ] ) 4 /it  WT 


(27) 


1 

TttWT 


between  impulses 


The  slight  differences  in  the  numerical  coefficients  can  be 

explained  but  we  will  not  digress  to  do  so  here.  The  interesting 

thing  to  observe  is  the  inverse  dependence  on  time-bandwidth 

product.  The  product  WT  is  proportional  to  the  number  of 

independent  samples  of  the  transmitted  process  x(t)  in  the  time 

interval  (-T,T) . Since  an  impulsive  impulse  response  does  not 

limit  the  bandwidth  of  the  reflected  signal  y(t),  the  product 

2 

x(t)y(t)  is  essentially  proportional  to  x (t)  and  still  has 

on  the  order  of  WT  independent  samples.  Equations  (7)  and  (9) 

indicate  that  Varj^  is  the  estimate  variance  with  the  noise 

n(t)  set  to  zero,  which  is  proportional  to  the  variance  of 

2 

essentially  an  average  of  a set  of  samples  of  x (t) . In 
normalized  form  this  varies  inversely  with  the  number  of 
independent  samples  involved,  i.e.,  WT.  The  result  can  be 
verified  by  approximating  the  integral  in  (9)  by  a sum  of 
independent  samples,  as  we  did  in  Section  II  B-l-b  above. 
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Now  compare  these  results  with  those  for  the  continuous 
impulse  response.  Equation  (19) . We  see  that  basically  where  W 
appears  in  the  impulsive  result,  appears  in  the  continuous 
result.  What  is  happening  is  this.  In  the  continuous  case, 
under  the  assumption  that  W the  bandwidth  of  the  reflected 

signal  y(t)  is  limited  to  W^  by  the  frequency  response  of  the 
target.  Thus  the  number  of  independent  samples  of  y(t)  available 
is  proportional  to  W^T,  not  WT,  and  this  change  carries  through 
to  the  result  (19) . Again,  this  inverse  W^T  dependence  can  be 
verified  by  a discrete  summation  approach,  although  it  is 
slightly  more  complicated  than  verifying  the  WT  dependence  in 
the  impulsive  case. 


b.  Var^  Term 

Sine  the  Varn  term  comes  about  from  cross-correlating 
n(t)  with  a replica  of  x(t),  it  does  not  depend  on  h(t) , and 
hence  in  un-normalized  form  is  unchanged  from  the  continuous 
case.  However,  the  normalized  form,  which  is  what  counts,  will 
change  in  accordance  with  (25),  and  becomes,  using  (22)  and  (25), 
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which,  using  (5),  can  be  written 


Var^[h(T) ] 


(E[fi(TQ) ]) 


ZI 

2 


E„/N 
^ n 


(29) 


where  Tq  represents  the  location  of  an  impulse  in  h(T), 

Comparing  (29)  with  (23)  we  see  that  the  growth  of 
variance  with  bandwidth  no  longer  occurs  in  the  impulsive  case. 
This  is  because  the  target  frequency  response  no  longer  limits 
the  bandwidth  of  the  reflected  signal  and  causes  energy  loss 
and  mismatch  at  the  receiver.  Again,  this  point  will  be  made 
more  explicitly  in  Section  II  C-2  below. 

C . Deterministic  Signal  Method 

The  basis  for  the  deterministic  signal  method  is 
shown  in  Figure  4.  The  transmitted  signal  s(t)  is  to  be 
thought  of  as  a real  "narrow"  pulse  with  energy  E given  by 


OD 

= / 


s ( t)  dt 


(30) 


To  remain  consistent  with  Section  B,  the  signal  will  be  assigned 
the  Gaussian  shape 


s(t)  = Sq  exp[-(2TTWt)  ] 


(31) 
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Fig,  4.  Deterministic  signal  method. 
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with  Fourier  transform 

S(f)  = Sq  exp[-(f/2W)^]  (32) 

where  one  can  show  that 

S ^ = 4t;W^S  ^ EW  (33) 

o o 

The  fact  that  the  signal  extends  to  infinity  and  that  the 
convolutions  to  follow  are  carried  out  with  infinite  limits  is 
not  significant.  The  same  analysis  has  been  carried  out  with 
signals  and  filters  having  finite  time  extent,  rectangular 
shapes  for  example,  with  no  change  in  the  basic  dependences, 
merely  form  factor  variations.  Truncating  the  Gaussian  shapes 
used  here  would  change  the  results  negligibly  while  introducing 
much  nuisance  value  into  the  analysis. 

The  reflected  signal  y(t)  is  s(t)  convolved  with  h(t), 
or 

s(C)h(t-C)dC  (34) 

00 

which  will  henceforth  be  denoted  by  s(t)  (^h(t).  The  received 
signal  is 

z(t)  = y(t)  + n(t)  (35) 
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which  is  then  passed  through  a receiver  filter  whose  impulse 
response  is  h„(t)  resulting  in  a waveform  which  acts  as  the 
estimate  of  h(t).  Thus, 


h(t)  = z(t)  0hj^(t)  = s(t)  0h(t)  0hj^{t)  + n(t)  0hj^(t) 


The  interference  n(t)  is  as  assumed  in  Section  B,  with  Gaussian 
shaped  spectrum  S^(f)  and  autocorrelation  Equations 

(1)  through  (3)  apply  with  Nj^  replacing  N^-  The  receiver  filter 
is  taken  to  be  a filter  matched  to  s(t),  and  therefore  has  an 
impulse  response  proportional  to  s(-t).  Thus  we  can  write 


hj^{t)  = hj^^  exp[-(27rWt) 


exp[-(f/2W) 


2 2 2 
^ 

Ro  Ro 


We  can  now  calculate  the  mean  and  variance  of  h(t)  for  the 
continuous  and  impulsive  impulse  response  models. 


1.  Continuous  Impulse  Response 


Assume  again  that  the  target  impulse  response  is 
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Gaussian  shaped  with  bandwidth  parameter  W^,  as  described  by 
Equations  (17)  and  (18).  From  (36)  we  can  write 

E[h(t)]  = s(t)  0h(t)  0h^(t)  = s(t)  0hj^(t)  @h(t)  (38) 

From  (31),  (32),  and  (37)  we  can  write 

s(t)  0h^(t)  = /Tn  exp  [- ( /2TTKt)  (39) 

For  large  W, 

s(t)©h„(t)  » So«Ro"o‘*^> 

Thus,  if  we  can  write  using  (38)  and  (40), 

E[h(t)  1 = SQHj^Qh(t)  (41) 

using  (36)  we  can  write 

oo 

Var[h(t)]  = E[m^(t)l  = = f ^^2) 

OO 

where  m(t)  is  the  noise  at  the  output  of  the  receiver  filter, 

R (t)  is  its  autocorrelation  and  S„(f)  is  its  power  spectrum, 
mm  m 

Therefore 
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Var[h(t)  ] 


(43) 


Sn(f)|H^(f) 


df 


This  evaluates  to 


Var  [h(t)  ] = /rf 


(44) 


and  using  (33)  and  (41),  letting  h^^l,  and  W>>Wj^,  we  can  write 


Var [h (o) ] 

/2  TT 

(E[h(o)  ])^ 

E/N„ 

To  understand  this  result,  consider  the  "narrowband" 
limit  W <<W^.  In  this  case  the  only  effect  of  convolving 
(s(t)  (^hj^(t))  with  h(t)  is  to  multiply  (s(t)  (2)hj^(t))  by  the 
amplitude  of  the  target  frequency  response  at  zero  frequency, 
namely  Thus  from  (38)  and  (39)  we  can  write 


E[h(o)]  = /27  (46) 

and  using  (44),  (46),  (33),  and  (18),  and  letting  h^=l,  we 

obtain  the  normalized  result 

(E(h(o)l)^  E/Nn 
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T 


I 

i 


I 

I 


Thus  the  ratio  of  the  normalized  variances  in  the  wideband  and 
narrowband  limits  is 


Var 

Var 


Wb 

nb 


(48) 


Two  effects  contribute  to  this  result.  In  the  narrowband  case 
the  receiver  is  matched  to  the  reflected  signal  y(t).  In  the 
wideband  case  the  bandwidth  of  the  reflected  signal  has  been 
limited  by  the  target,  while  the  receiver  operates  wideband  to 
match  the  transmitted  signal.  Thus  excess  noise  enters  the 
receiver  relative  to  the  matched  case,  the  effect  being 
proportional  to  W/W^.  The  other  effect  is  that  the  target 
rejects  most  of  the  transmitted  energy  outside  its  nominal 
bandwidth  Wj^.  Energy  must  be  transmitted  "out  of  band"  in 
order  for  the  expected  pattern  of  the  estimate  to  contain  the 
desired  detail  or  resolution,  but  most  of  this  energy  is  lost 
from  the  standpoint  of  the  variability  of  the  estimate.  This 
effect  is  also  proportional  to  W/Wj^. 

We  see  then  that  "too  much"  bandwidth  is  a bad  thing, 
using  either  the  deterministic  or  the  random  signal.  If  you 
ask  for  more  resolution  or  detail  than  you  really  need,  a 
noisier  estimate  results. 
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2.  Impulsive  Impulse  Response 


Assume  again  that  the  target  impulse  response  is  not 
a continuous  curve,  but  contains  impulses.  For  simplicity,  the 
single  impulse  or  point  target  model  will  be  treated.  Thus 
h(t)  is  taken  to  be 


h(t)  = h^u^(t) 


Then  from  (38)  we  can  write 


E[h(t)]  = h^s(t)  0hj^(t) 


(49) 


(50) 


which  says  that  the  expected  shape  of  the  estimate,  ideally  an 
impulse,  is  the  shape  of  the  pulse  at  the  output  of  the  receiver 
filter.  Applying  (39)  and  (33)  we  get 


(E[h(o)])2  = 2tt  = /?¥  EWHj^^^h^^  (51) 


The  un-normalized  variance  is  still  given  by  (44) , and  the 
resulting  normalized  expression  is 

Var [h (o) ] _ ^ 

(E[h(o)])^  2 


(52) 

E/N 

n 


Ill . 


Conclusions 


1 


In  Section  I we  addressed  heuristically  the  question 
of  whether  the  use  of  random  rather  than  deterministic  signals 
had  any  advantage  for  impulse  response  estimation.  In  Section 
II  we  defined  an  objective  basis  for  making  the  comparison, 
carried  the  comparison  out,  and  developed  some  insight  into  the 
results.  In  this  section  we  will  summarize  what  has  been 
learned  and  briefly  comment  on  it. 

The  basic  performance  expressions  derived  are  brought 
together  in  Table  1 below. 


TABLE  1 

NORMALIZED  VARIANCE  OF  IMPULSE  RESPONSE 
ESTIMATE 


RANDOM 

SIGNAL 


CONTINUOUS  IMPULSE  RESPONSE 


VARIANCE  DUE 
TO  SIGNAL 


1 1 
2 Wj^T 


DETERMINISTIC] 
SIGNAL  ' 


NONE 


VARIANCE  DUE 
TO 

INTERFERENCE 


/2tt 


W 

E/N 


/2i 


W 

E/N 


IMPULSIVE  IMPULSE  RESPONSE 


VARIANCE  DUE! 
TO  SIGNAL 


1 

277  WT 


NONE 


VARIANCE  DUE 
TO 

INTERFERENCE 


/T  1 

2 E/N 


n 


1 

2 E/N, 
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These  results  apply  in  the  limits  on  time  and  bandwidth  pre- 
viously discussed,  and  the  numerical  coefficients  correspond 
to  the  Gaussian  shape  assumptions  and  definition  of  bandwidth 
made.  The  important  point  is  that  bandwidth  has  been  defined 
consistently  in  the  sense  that  a given  value  of  W implies 
exactly  the  same  resolution  or  detail  in  the  expected  estimate 
in  either  the  deterministic  or  the  random  case. 

The  table  shows  clearly  what  the  answer  is  regarding 
the  central  question  addressed  in  this  report,  namely  the 
relative  ability  of  the  two  methods  to  suppress  uncorrelated 
interference.  Their  performance  is  identical  for  either  type 
of  impulse  response  model  when  they  employ  equal  resources  such 
as  bandwidth,  average  energy  consumption,  etc. 

As  described,  the  random  signal  method  experiences  an 
additional  variance  component,  although  this  can  be  reduced  to 
a desired  level  by  adjusting  the  parameters.  Alternately  one 
can  select  a particular  realization  of  the  random  process  and 
use  it  repeatedly  as  a radar  signal  just  like  any  other 
deterministic  signal,  e.g.,  LFM,  binary  phase  coded,  etc. 

Finally,  we  have  pointed  out  that  with  either  type  of 
signal  the  use  of  "excessive"  bandwidth  carries  with  it  the 
penalty  of  a more  variable  estimate.  This  should  be  kept  in 
mind  in  discrimination  applications  where  achieving  fine  resolution 
on  the  shortest  targets  in  question  may  in  fact  be  unnecessary  or 
even  undesirable* 
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APPENDIX 

The  purpose  of  the  appendix  is  to  describe  in  somewhat 
more  detail  than  was  appropriate  in  Section  II  how  some  of  the 
variance  expressions  are  arrived  at.  In  the  random  signal  case 
we  start  with  equation  (14)  giving  the  definition  of  the 
estimator  variance  as 


Var[h(T)  ] = E[(h(T)  - h(T)  ) 


Using  the  definition  of  h(T)  from  Equation  (9), 


h (t)  = 


T 

^ / 

TN 

X -T 


X ( t) z ( t+  t)  dt 


the  expression  given  in  (11)  for  the  expected  value  of  h(T), 


E[h(T)]  = 


1 T 

v=-T-^ 


( v)  h ( v+t)  dv  dC 


and  the  signal  relationships  from  (6)  and  (7) , 


z(t)  = y(t)  + n(t) 


and  y(t)  = 


x(C)  h(t-5)  d^ 


we  can  write 


T T-  f 

h(i)  - fi{T)  = 2TNx  f f 

C=-T  v=-T-C 


[x(5)x(5+v)  - ( v)  ] h ( v+t)  dvd^ 


1’ 

^x  f 


X ( p)  n ( p+t)  dp 


Think  of  the  first  (double)  integral  as  a term  "A"  and  the 

second  integral  as  a term  "B".  Then  the  square  of  the  left 

2 2 2 

side  of  (A6)  is  equal  to  (A+B)  which  equals  A + 2AB  + B . 
2 

The  "A  " term  is  the  square  of  the  double  integral  and  can 
therefore  be  written 


Ty.  T-r  T T-t 


a2  = 


4T^N 


f f f f 


[x(C)x(C+v)-Rxj^(v)  ] 


X 5=-T  v=-T-5  t=-T  p=-T-t 


[x(t)  x(t+p) -R^^(p)  1 h(v+T)h(p+T)  d^dvdtdp 


2 

The  "B  " term  is  the  square  of  the  single  integral  and  can 
therefore  be  written 


b2  = 


-T 

Y — 2 j/  X ( p)  x(  t)  n ( p+x)  n(  t+x)  dpdt 

N 'U'  m 


42 


The  cross  term  "2AB"  is  the  product  of  the  single  and  double 


t 


integrals,  which  can  be  written 


2AB  = 


2T‘ 


T T-f  T 

^ / / / 


[x(^)x(r  + v)-R  ^(v)  ] 


^=-T  v=-T-^  p=-T 


(A9; 


X ( p) n ( p+T ) h { v+T ) dpdvd^ 


The  expressions  (A7)  through  (A*))  comprise  the  quantity 
[h(T)  - E(h{T)]'^.  The  variance  of  h{T)  is  then  the  sum  of 
their  expected  values.  In  each  case,  the  expectation  is  moved 
inside  the  integral.  The  first  thing  we  note  is  that  because  of 
the  independence  of  x and  n,  the  entire  integrand  in  (A9)  is 
proportional  to  E[n(p+T)]  which  is  zero  (see  Equation  (8),  Section 
II  B) . Thus  (A9)  contributes  nothing,  and  we  have  only  (A7) 
and  (A8)  to  consider. 

The  expected  value  of  (A8)  is  the  Var^  term.  If  we  make 
the  change  of  variables  C=  p-t,  we  can  rewrite  (A8)  as 


Var^ [h (t) ] 


1 

2 2 


^=-T-t 


E[x(t)x(t+5)n(t+T)n(t+T+f;)  ] 

(AlO) 


dCdt 
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Again  using  (8)  because  x and  n are  independent,  we  arrive  at 


Varj^[h(T)  ] 


4T^N 


X 


t=-T  C=-T-t 


(All) 


This  is  the  basic  result  for  Var  . If  we  substitute  the 

n 

particular  Gaussian  shapes  assumed  in  Section  II  for  and 

Rj^^,  the  integration  can  be  carried  out  in  closed  form,  resulting 
in  an  expression  involving  the  error  function  "erf."  When 
this  expression  is  evaluated  for  large  WT,  expression  (22)  for 
Var^^  results. 

The  expected  value  of  (A7)  is  the  Var^  term.  The  pro- 
cedure is  as  follows.  The  expectation  is  moved  inside  the  inte- 
gral and  the  integrand  is  multiplied  out.  Three  of  the  four 
terms  that  result  are  of  the  form  R (v)R  (p ) h (v+x ) h (p+x ) . The 

XX  XX 

fourth  term  is  of  the  form  E [x  (5)  x (^-t-v)  x (t)  x (t+p)]  h (v+x ) h (p+x ) . 

The  expectation  can  be  calculated  by  applying  the  "moment  theorem" 
for  Gaussian  variables",  which  can  be  stated  in  the  form 


E[a6Yi5]  = E[a6]  E[y6]  + E[aY]  E[65]  + E[a6]  E [ By] 

(A12) 
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where  a,B,Yf  and  6 are  real  zero  mean  Gaussian  random  variables. 
If  {A12)  is  applied,  one  of  the  three  terms  that  results  again 
contains  the  factor  (v ) Rxx  (P ) / and  just  cancels  the  previous 
terms  containing  that  factor.  Thus  all  that  remains  are  the 
other  two  terms  that  result  from  applying  (A12) , and  we  arrive 
at 


Varx[h(T)]  = 


4T^N 


X 


m rn_  m 

/ /'  / 

; =-T  v)=-T-r  1-  = -' 


[Rj^^(C-t)Rj^^(C-t+v-p)  + Rj^jj(t+p-C)Rxx(-+^"t)l 


(A13) 


h (v+t) h (p+t) dpdtdvdC 


This  is  the  basic  result  for  Varxf  and  it  is  not  very  informative 
as  it  stands.  It  remains  to  turn  it  into  a more  useful  form 
in  both  the  continuous  and  impulsive  impulse  response  cases. 

In  the  continuous  case  the  assumption  that  R is 

XX 

narrow  relative  to  h,(WTj^  >>  1),  allows  us  to  make  the  substitution 
R^x  ( ’ ) ^ ^x'^o  ^ ^ (A13)  . If  we  do  this  and  then  integrate 

first  over  p and  then  over  t,  the  result  is 


A 

Varx [h (t) ] 


7 

E=- 


7 

v=-T-5 


h (v+t) [h (v+t) 


+ h (T-v) ] 


dvdv 

(A14) 
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and  if  we  assume  for  convenience  that  h is  even  and  evaluate 
at  the  typical  point  t=0,  the  result  becomes 


T.  T-C 


Var^[h(o)]  = 


(v)dv  (A15) 


^=-T  v=-T-^ 


which  is  Equation  (16)  of  Section  II. 

In  the  impulsive  case,  we  have  by  assumption  that  h is 
impulsive  relative  to  Rxx'  whereas  in  the  continuous  case  we 
treated  R ^ as  impulsive  relative  to  h in  the  limit  of  interest. 
If  we  return  to  (A13) , integrate  first  on  v,  then  on  p,  and 
assume  here  for  simplicity  the  point  target  model  h(T)  = 
hoUo(T),  we  obtain 


Var  [h(T)]  = — 2 2 II 


dCdt  (A16) 


If  we  evaluate  at  the  typical  point  t=0,  and  change  variables 


we  arrive  at 


J J 

^ 5=-T  n=-T-C 


Rj^x  dndC  (A17) 


This  can  be  evaluated  just  like  (All),  and  Equation  (26)  results, 
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